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ABSTRACT 

Aims. Massive stars are known to demonstrate significant spectroscopic and photometric variability over a wide range of timescales. 
However the physical mechanisms driving this behaviour remain poorly understood. Westerlund 1 presents an ideal laboratory for 
studying these processes in a rich, coeval population of post-main sequence stars and we present a pathfinding study aimed at charac- 
terising their variability. 

Methods. To do this we utilised the large body of spectroscopic and photometric data that has accumulated for Wdl during the past 
decade of intensive studies, supplemented with the sparser historical observations extending back to the early 1960s. 
Results. Despite the heterogeneous nature of this dataset, we were able to identify both spectroscopic and photometric variability 
amongst every class of evolved massive star present within Wd 1. Spectroscopic variability attributable to both wind asphericity and 
photospheric pulsations was present amongst both the hot and cool hypergiants and the former, also with the Wolf Rayets. Given the 
limitations imposed by the data, we were unable to determine the physical origin of the wind structure inferred for the OB supergiants, 
noting that it was present in both single pulsating and binary stars. In contrast we suspect that the inhomogineities in the winds of the 
Wolf Rayets are driven by binary interactions and, conversely, by pulsations in at least one of the cool hypergiants. Photospheric pul- 
sations were found for stars ranging from spectral types as early as 09 I through to the mid F Ia + yellow hypergiants - with a possible 
dependence on the luminosity class amongst the hot supergiants. The spectroscopically variable red supergiants (M2-5 la) are also 
potential pulsators but require further observations to confirm this hypothesis. Given these findings it was therefore rather surprising 
that, with the exception of W243, no evidence of the characteristic excursions of both luminous blue variables and yellow hypergiants 
was found. Nevertheless, future determination of the amplitude and periodicity of these pulsations as a function of temperature, lu- 
minosity and evolutionary state holds out the tantalising possibility of constraining the nature of the physical mechanisms driving 
the instabilities that constrain and define stellar evolution in the upper reaches of the HR diagram. Relating to this, the lack of secu- 
lar evolution amongst the cool hypergiants and the presence of both high-luminosity yellow hypergiants and red supergiants within 
Wdl potentially place strong constraints on post-main sequence evolutionary pathways, with the latter result apparently contradicting 
current theoretical predictions for > 25M G stars at solar metallicites. 
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C"" i 1. Introduction stars, and instead that the H-rich mantle is shed in short-lived im- 
pulsive events associated with the 'zoo' of transitional objects. 

O . A large body of observational evidence exists that massive post- Inde ed observational evidence of this hypothesis is provided by 

main sequence (MS) stars are subject to a variety of (pulsational) the p reS ence of fossil ejecta surrounding luminous blue vari- 

instabihties (e.g. Humphreys & Davidson 1994 Sterken [1977J ables (LBVs : e.g. Clark et al. l2003l . yellow hypergiants (YHGs; 

. r-j . Burki [19781 de Jager [19981 Wood et al. [1983). Indeed, it has C astro-Carrizo et al. l2007l l and red supergiants (RSGs; Schuster 

[XJ ■ been supposed that such instabilities contribute to the absence of et al |2QQ6l > and the direct association of dramatic increases in 

H ; cool (T< lOkK) hypergiants above log(L/L )~5.8, giving rise to mass loss with e pi so des of instability in both hot and cool lumi- 

the the Humphreys Davidson limit that spans the upper reaches nous stars Carand p Cas, respectively; Smith & Gehrz lT998l 

of the HR diagram. Lo b e i et al [20Tj3l >. Thus, determining the duration, duty cycle 

Characterisation of such instabilities is important for two an d mass i oss rate f such events appears to be necessary for the 

reasons. Firstly, one might hope to use the pulsational properties construction of a theory of massive stellar evolution, 
to investigate the internal structure of massive stars via astero- 

seismology. Secondly, there is increasing reason to suppose that These instabilities are manifest in both photometric and 

MS mass loss is insufficient to permit the transition of high-mass spectroscopic variability, but other physical processes, such as 

stars from unevolved O-type objects to H-depleted Wolf Rayet large and small scale wind inhomogineities, also lead to simi- 
lar behaviour. Line driven winds are intrinsically unstable, lead- 
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* Based on observations collected at the European Southern namical simulatons by Cranmer & Owocki H1996JI demonstrate 

Observatory, Paranal and La Silla Observatories under programme IDs that large-scale photospheric perturbations - originating in non- 

ESO 067.D-0211, 069.D-0039, 071.D-0151, 073.D-0327, 075.D-0388, radial pulsations and/or magnetic fields - can lead to large scale 

081.D-0324 and 383.D-0633. rotating wind structures. Moreover, spectrophotometric variabil- 
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ity is an important observational signature of wind collision 
zones in massive binaries (e.g. Lewis et al. 119931 Stevens & 
Howarth 1999 ). Of these phenomena, the degree of wind clump- 
ing present in massive stellar winds is critical to the observa- 
tional determination of the mass loss rate, while the characteri- 
sation of the properties of the massive binary population is im- 
portant for fields as diverse as star formation, the nature of SNe 
progenitors and the formation channels of X-ray binaries. 

Consequently, considerable observational effort has been ex- 
pended in the identification and characterisation of variability 
in massive (binary) stars. Such studies have typically centred 
on field stars, with the attendant difficulties of the determina- 
tion of stellar luminosity and hence mass and age. However the 
development of multiplexing spectrographs and the identifica- 
tion of a large number of massive young Galactic clusters - of 
which Westerlund 1 is an examplar (henceforth Wdl) - presents 
a unique opportunity to advance this field, offering a large co- 
eval population of massive stars of identical initial metallicity 
and well constrained masses. 

First identified by Westerlund d 1 96 It . Wdl was subject to 
sporadic optical and near-IR photometric studies in the follow- 
ing four decades (Sect. 2), with the high reddening towards it 
preventing a spectroscopic study until that of Westerlund ( 1987; 
West87). However, despite the discovery of a remarkable popu- 
lation of both early and late super/hypergiants, the cluster lan- 
guished in relative obscurity until the turn of the century, when 
modern spectroscopic studies - motivated by the unexpected ra- 
dio properties of Wd 1 - first identified a rich population of Wolf- 
Rayets (WRs; Clark & Negueruela 2002) and subsequently re- 
vealed Wdl to be the first example of a super star cluster iden- 
tified in the Galaxy (Clark et al. l2005al henceforth C05). In the 
following decade, Wdl has been the subject of a large number 
of studies from X-ray to radio wavelengths aimed at determin- 
ing both bulk cluster properties and those of the individual con- 
stituent stars. The combined dataset therefore offers the potential 
of investigating the variability of cluster members over a signif- 
icant time period (>40 yr for the brightest stars; Sect. 2) and 
we present the results of such an analysis in this manuscript. In 
Sect. 2 we present and describe the reduction of the new data 
included in this work and discuss the utilisation of data available 
in previous studies. In Sect. 3-7 we analyse the available data 
as a function of evolutionary state and provide a discussion and 
summary in Sect. 8. 

2. Observations & Data Reduction 

Between its discovery in 1961 and the era of intensive observa- 
tions post 2001, Wd 1 was the subject of a number of photo- 
metric and spectroscopic studies ranging from optical to mid-IR 
wavelengths; these are summarised in Table 1. Unfortunately, in 
many cases observations collected in single works were not con- 
temporaneous, complicating analysis; an important example is 
the photometry and spectroscopy presented by West87, which 
were obtained 15 years apart. An additional difficulty in com- 
paring multiepoch photometry is the different filter and detector 
(photographic plate, photometer and CCD) combinations em- 
ployed by the various studies, for which accurate spectral re- 
sponses were unavailable. Combined with different extraction 
techniques (important due to the crowded nature of the cluster) 
and the lack of appropiate red standards, we found the number of 
observationally robust results that could be drawn from the his- 
torical data were limited. Finally, the recent near-IR photometric 
observations of Brandner et al. (2008 ) provide no extra informa- 
tion due to the saturation of the near-IR bright supergiants. 
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Fig. 1. Comparison of the photometric datasets presented by 
C05, Bonanos (120071) and West87. Note the zero point offsets 
present between all three datasets. In all three panels outliers are 
flagged, along with observations of such stars in the other wave- 
bands, where available. Also note that, for a given star, photom- 
etry is not always available for all three datsets (Sect. 2.1). 
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w# 


1966 


1973 


Spectral Classification 
1973-4? 1981 1995 


2002 


2006 


4 


OB 


OB 




GO la 


OB 


F2Ia 


<G5 la 


6 


OB 


- 




Bl V 




B0.5 lab 


- 


7 


OB 


OB 




AO la 


OB 


B5 la 


- 


8a 


OB 


- 




GO la 


OB 


F5 la 


<G5 la 


9 


OB 


7 


? 


Be 


OB 


sgB[e] 


<G5 la 


10 


OB 


- 




BO III 




B0.5I+OB 


- 


11 


OB 


- 




Bl II 




B2 la 


- 


12a 


OB 


OB 




A2 la 


OB 


A5 la 


<G5Ia 


16a 


OB 


OB 




A2 la 


OB 


A2 la 


- 


20 


M 


M 


M 


M6 la 




<M6I 


M5 la 


26 


M 


M 


M 


M2 la 


M 


<M6I 


M5-6Ia 


28 


OB 


- 




BO II 




B2 la 


- 


31 


OB 


- 




Be 




B0I+OB 


- 


32 


OB 


OB 




G5 la 


OB 


F5 la 


<G5 la 


33 


OB 






B8Iab 


OB 


B5 la 




70 


OB 






B8Iab 




B3 la 




71 


OB 






B8Iab 




B2.5 la 




75 


M 


M 


M 








M4Ia 


78 


OB 






Bl la 




Bl la 




237 


M 


M 


M 


M6 




<M6I 


M3 la 


243 


OB 


OB 




B2Ia 




A2 la 




265 


OB 


OB 




GO la 




F5 la 


<G5 la 



Note that classifications given in itallics have been derived from photometric data alone, while the 2006 classifications are derived from 2.3/im 
spectroscopy (Mengel & Tacconi-Garman 2007 , 2008 ). Crowther et al. ( 2006 ) provide a comparable table for the WR population. 



Table 3. Summary of the properties of the spectroscopic datasets employed in this study 



Date 


MJD 


Telescope 


Instrument 


Setting 


A 

(A) 


Resin. 


Reference 


23-25/06/01 


52083-5 


ESO 1.52m 


Boiler & Chivens 


Grat. #1 


6000-11000 


500 


C05 


07/06/02 + 
06-08/06/03 


52432 
52796-8 


NTT 


EMMI 


Grat. #7 
Grat.#6 


6310-7835 
8225-8900 


2600 
5000 


C05 

Negueruela & Clark ( 2005 1 


12-13/06/04 


53168-69 


VLT 


FORS longslit 

FORS MXU 
FORS MXU 


G1028z 
G1200R 
G1200R 
G1028z 


7730-9480 
5750-7310 
7730-9480" 
5750-7310" 


-7000 
-7000 
2560 
2140 


Negueruela et al. dsubmittedb 


25/03/05 + 
29/05/05 + 
13/07/05 


53454 
53519 
53564 


VLT 


FLAMES + 
MEDUSA 


LR6 
LR8 


6438-7184 
8206-9400 


8600 
6500 


Ritchie et al. d2009bt 


29-30/06/05 


53550-1 


NTT 


SOFI 


IJ Grism 
HK Grism 


9500-13500 
15000-22700 


1000 
1000 


Crowther et al. (|2006|» 


17/02/06 


53783 


NTT 


EMMI 


Grat. #6 


5740-8730 


1500 


Marco & Negueruela ( 2009 ) 


11-12/03/06 


53805-6 


VLT 


ISAAC 


SWS(MR) 


22490-23730 


9000 


Mengel & 

Tacconi-Garman ( 2007 ) 


20+29/06/08 + 
18+24/07/08 + 


54367+76 
54665+71 


VLT 


FLAMES + 
MEDUSA 


HR21 


8484-9001 


16200 


Ritchie et al. d2009a& 



14+17/08/08 + 54692+95 

04+15/09/08 + 54713+24 

19+25/09/08 + 54728+34 

14+18/05/09 + 54965+69 
20/08/09 55063 



"Note that the nominal wavelength range is given for the 2004 FORS MXU observations; in practice this can be displaced by a few hundred A 
depending on the position of an individual star on the CCD. 
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Table 4. Summary of spectroscopic data for supergiants of spectral type B2 and earlier. 



ID 


Spectral 


RA, DEC (J2000) 


R 


I 


1.5m 


NTT 


NTT 


VLT 


VLT 


Notes 




Type 








(2001) 


(2002) 


(2003) 


(2004) 


(2008/9) 




W2a 


B2Ia 


16 46 59.7 -45 50 51.1 


14.23 


11.73 


• 


#6 




R+z 




A, cand. binary, Ho- var. 


W6a 


B0.5Iab 


16 47 04.0 -45 50 21.0 


15.80 


13.16 


• 




#7 


R+z 


• 


P, cand. binary, Ho- var. 


W8b 


B1.5Ia 


16 47 06.0 -45 50 25.0 


- 


- 








R+z 


• 


pulsator 


W10 


B0.5+OB 


16 47 03.3 -45 50 34.7 


- 


- 


• 






R+z 




A, binary 


Wll 


B2Ia 


16 47 02.2 -45 50 47.0 


14.52 


11.91 




#6+#7 




R+z 




A 


W15 


091b 


16 47 07.0 -45 50 28.0 


16.38 


13.75 


• 






R 


• 


A 


W19 


Blla 


16 47 04.9 -45 50 59.1 


15.21 


12.37 




#6+#7 




R+z 




A, Ha var. 


W21 


B0.5Ia 


16 47 02.0 -45 51 11.0 


15.56 


12.74 








R+z 


• 


A, pulsator 


W23a 


B2Ia+OB 


16 47 03.0 -45 51 08.0 


14.91 


12.07 




#6+#7 




R+z 


• 


A, binary, Ho- var., pulsator 


W24 


09Iab 


16 47 03.0 -45 51 11.0 


15.96 


13.24 




#6+#7 




z 


• 


cand. pulsator 


W28 


B2Ia 


16 47 04.7 -45 50 38.4 


14.26 


11.64 




#6+#7 




R+z 




A, Ha var. 


W29 


091b 


16 47 04.4 -45 50 39.8 


16.02 


13.38 


• 






R+z 






W30a 


O+O 


16 47 05.0 -45 50 39.0 


15.80 


13.20 




#6+#7 




R+z 


• 


binary, Ha var. 


W35 


09Iab 


16 47 04.2 -45 50 53.5 


16.00 


13.31 


• 






R+z 






WA 1 


(J91ab 


16 4/ U2. / -43 30 36.9 


13.39 


1 O "70 

12. la 


• 






R+z 






W43a 


Blla 


16 46 03.5 -45 50 57.3 


15.22 


12.26 


• 


#6+#7 




R+z 


• 


A, binary, Ha var. 


W55 


BOIa 


16 46 59.0 -45 51 29.0 


15.25 


12.67 


• 








• 


A 


W60 


BOIab 


16 47 05.0 -45 51 51.0 


15.96 


13.28 




#6+#7 




R+z 


• 


A 


W61a 


B0.5Ia 


16 47 02.3 -45 51 41.6 


14.62 


12.01 




#6+#7 




R+z 




A, Ha var. 


W61b 


09.5Iab 


16 47 02.6 -45 51 41.6 


16.00 


13.31 




#6+#7 




z 






W74 


09.5Iab 


16 47 07.0 -45 50 10.0 












R+z 


• 


A 


W78 


Blla 


16 47 02.0 -45 49 56.0 


14.54 


12.04 








R+z 


• 


A, pulsator 


W238 


Bllab 


16 47 05.0 -45 52 27.0 


14.98 


12.45 


• 






R+z 


• 


A 



Note that W2a was also observed in 2005 March and W61a in 2005 May+July (VLT FLAMES+Medusa). Spectral classifications are from 
Negueruela et al. dsubmittedl l. The final column indicates whether a given star is identified as a Periodic or Aperiodic photometric variable, 
(candidate) binary, Ha variable and/or (candidate) photospheric pulsator (Sect. 3). 



Table 1. Summary of published observations utilised in this 
study 



Year 


Dataset 


Reference 


1961 


VI band phot. 


Borgman et al. ( 1970) 


1966 


VRI band phot. 


West87 


1969 


K band phot. 


Borgman et al. ( 1970 1 


1973 


V+narrow band phot. 


Lockwood dl974b 


1973-4 


2-20yum phot. 


Koornneef (fT977l 


1995 


VI band phot. 


Piatti et al. (fl998l 


2001 


BVRI band phot 


C05 


2006 


BVRI band phot. 


Bonanos ( 2007 ) 


1981 


RI band spec. 


West87 


2001 


RI band spec. 


C05 


2002 


RI band spec. 


C05 


2003 


RI band spec. 


Negueruela 






& Clark d2005b 


2005 


IJ & HK spec. 


Crowther et al. (|2006) 


2006 


K band spec. 


Mengel & 






Tacconi-Garman (|2007) 



2.1. Photometric data 

Comparative analysis of multiepoch photometric data obtained 
prior to 2001 was not attempted for the above reasons and, with 
the exception of West87, the limited number of stars observed 
in such works, which precluded the determination of empirical 
offsets between different datasets (e.g. Fig.lfl Hence the only 
robust results to be drawn from these observations are the coarse 
spectral classifications derived from colour indices for the (con- 

1 Note that the limited stellar identifications and lack of co-ordinates 
presented by Piatti et al. ( 1998 1 limit the utility of those observations. 



temporaneous) datasets presented by West87, Lockwood ( 1974 ), 
Koornneef d 19771 1 and Piatti et al. (119981 1. Unfortunately these 
just permit the classification of stars into OB and M supergiant 
classifications, for which we found the former to contain stars 
of spectral type from late O to mid GIa + . These results are sum- 
marised in Table 2. 

These problems persist for modern, post-2001 datasets. This 
is apparent in the systematic offset between the V band datasets 
present in C05 and Bonanos (2007; Fig. 1) due to the difficulty 
in determining the zero point offset because of the lack of ap- 
propiate red standards^. However the large number of stars in 
common between West87, C05 and Bonanos (120071 ) enabled the 
offsets to be determined empirically and hence a direct compar- 
ison to be made between the three epochs (Fig. 1). 

This approach allows potentially variable outliers to be iden- 
tified, albeit with three corrolaries. Firstly observations of a par- 
ticular star are not always present in each dataset; for example 
W4, 20, 26, 32, 42, 75 & 243 are absent from C05. Secondly, 
observations of individual objects are likewise not present in all 
3 wavebands of a given dataset. For example, comparable pho- 
tometry for W243 between West 87 and Bonanos (2007 ) is only 
available in the V band, thus the apparent discrepancy between 
the two epochs is impossible to verifyQ (complicated by the fact 
that the empirical offset is ill defined for stars of this magni- 
tude). Finally, issues arising from filter and photographic emul- 

2 Bonanos (2007) successfully made such a correction for R and I 
band, but the methodology employed - comparison of the brightest 100 
stars in common between the 2 datasets - was not completely successful 
for the V band data. We suspect this to be because the brightest stars 
in the V band were dominated by foreground stars of relatively low 
reddening in comparison to the stars within Wdl. 

3 The same is the case for W44 & 57a, and W57a in the comparison 
between West87 and C05 datasets. 
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sion response to very red objects (West87, Bessell 1986) cast 
doubt upon the reality of the offsets between the data presented 
by West87 for the RSGs W26 & 237 and other epochs, although 
this would not explain the apparent variability of W237 in V & 
R bands between 2001-6 which we consider likely to be of as- 
trophysical origin (Sect. 7). 

Given the above, in addition to W237 we may only tenta- 
tively identify a handful of stars - W9, 57a & 243 - as poten- 
tially variable, although with only 3 epochs of observations over 
a 20 yr baseline and a lower limit on detectable variations of 
>0.5mag (based on the scatter between datasets evident in Fig. 
1), these are hardly stringent constraints. Nevertheless, foreshad- 
owing the detailed results and discussion presented in Sects. 
3-7, both hot and cool transitional stars may demonstrate long 
term (>decades) variability of many magnitudes amplitude; we 
would have expected to identify such behaviour photometrically 
if present within the evolved stellar population of Wdl . 

2.2. Spectroscopic data 

Prior to the discovery of a rich population of WRs within Wdl, 
the only spectroscopic survey to have been undertaken was that 
of West87. Motivated by this discovery, several optical and near- 
IR studies have subsequently been made (Table 1); the published 
results of all these works are summarised in Table 2. Given the 
varied science goals of the proposals from which these obser- 
vations have been drawn, the resulting spectra sample a num- 
ber of different wavelength ranges and corresponding resolu- 
tions; these are summarised in Table 3 with, where available, 
a reference for a more detailed description including the reduc- 
tion procedure adopted. Of these data, those obtained in 2003 
June (NTT/EMMI), 2005 March, May & July (VLT/FLAMES), 
2006 February (NTT/EMMI) and 2008 September & 2009 May 
& August (VLT/FLAMES) are published here for the first 
time, as are the observations of the YHGs from 2004 June 
( VLT/FORS) with the corresponding data for the OB stars found 
in Negueruela et al. (submitted). In total this amounts to in ex- 
cess of a hundred new spectra. 

Special mention must be made of the optical spectra obtained 
in 2001 June (C05) and the near-IR spectra obtained in 2006 
March (Mengel & Tacconi-Garman I2U071 120081 . The former 
are of significantly lower resolution and S/N than later epochs. 
Following the analysis presented in Negueruela et al. (submitted) 
we may utilise these data to adequately distinguish between ex- 
treme B hypergiants/WNL stars (Group 1 stars following C05, 
with the additional criteria that Ha must be in emission), late 
O/early B supergiants (Group 2; 09-B2 I), mid-late B super- 
/hypergiants (Group 3; B2.5-9I (+) ) and YHGs. Thus, while these 
data are inadequate for the analysis of line profile variability 
they do permit the gross spectral classification of cluster mem- 
bers, hence extending the baseline for the identification of e.g. 
the characteristic temperature excursions of LBVs by a further 
year. 

The latter spectra sample a relatively poorly studied region of 
the near-IR waveband, for which a comprehensive spectral atlas 
has yet to be constructed. Nevertheless, the relative strength of 
the CO bandheads may permit classification of cool supergiants 
(Figer et al. 120061 Yamamuro et al. 120071 1. Mengel & Tacconi- 
Garman (2007 2008) use this diagnostic to classify the 4 RSGs 
within Wdl, while the lack of CO bandhead absorption features 



4 However they are sufficient to identify dramatic changes in line pro- 
files such as the appearance of He emission in W23a and W43a; C05 
and Sect. 3. 



in the spectra of W4, 8a, 9, 12, 32 & 265 suggest a spectral type 
of G5 la or earlier. 

Summaries of the data available for late O/early B super- 
giants and the more evolved, transitional stars are presented in 
Tables 4 & 5 respectively. A comparable table for the WRs may 
be found in Crowther et al. (2006), noting that the three WRs in- 
cluded in Table 5 are the sole objects for which multiple epochs 
of high resolution data are available, thus permitting analysis of 
line profile evolution. Given the recent comprehensive analysis 
by Ritchie et al. d2009bl) . the LBV W243 has been omitted from 
this study. 



Table 5. Summary of the spectroscopic datasets available for B2 and later super-/hypergiants and Wolf-Rayets. 



ID 


Spect. Type 


RA, DEC (J2000) 


R 


I 


06/01 06/02 


06/03 


06/04 


03/05 


05+07/05 


02/06 


03/06 


2008/9 


Notes 


W4 


F3 Ia + 


16 47 01 42 -4*5 50 37 1 


11.18 


9.15 


. #6+#7 


#6 




LR6+8 


LR6+8 








A,P+R 


W5 (S) 


WNT /RTa + 

VV 1 > 1^ IJ I <-l 


1 6 46 07 97 -45 "SO 1 9 5 


14.98 


12.48 




#7 


R+z 


LR8 












W7 


B5 Ia + 


16 46 03 62 -45 50 14 2 


12.73 


9.99 


a #6+#7 


#6 +#7 






LR6+8 








A 


W8a 


F8 Ia + 


1 6 47 04 79 -45 50 24 9 


12.64 


9.89 


a #6+#7 








LR6+8 










W9 


orrRrpl 

sg D L C J 
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Note the spectral classification of the OB stars are from Negueruela et al. (submitted I, the YHGs from this work (Sect. 6) and the RSGs from Mengel & Tacconi-Garman (2007 1. Where appropriate, 
the Notes column lists stars as Eclipsing or Aperiodic photometric variables (Bonanos 20071, X-ray (Clark et al. 2008 1 and Point and/or Resolved radio (Dougherty et al. l2010t sources. Finally, the 
binary status of 5 stars - as determined by Ritchie et al. ( 2009a ) from multiepoch radial velocity spectroscopic data - is presented, noting that insufficent data currently exist to confirm the status of 
the other stars (but note Sect. 6 for W265). 
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Fig. 2. Left panel: Ha for W2a (B2 la) from 2004 - the line is 
strongly in emission in this spectrum but was absent in the low 
S/N spectrum from 2001 (not shown). Right panel: Ha profile 
for W6a (B0. 5 lab) in 2003-4; note the temperature dependent 
He lAA 6678, 7065 A lines also appear to vary in strength. 
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Fig. 4. LPV in the 2008-9 VLT/FLAMES spectra Cm+Pal6 
blend in W24; in comparison little variability is observed in the 
unblended Pall profile. Note that the strength of Cm line in 
the 2002 NTT and 2004 VLT spectrum, which are not shown 
for clarity, are comparable to the strong (e.g. MJD 54665.1) and 
weak (e.g. MJD 54637.2) 'states' respectively. 



3. The 09-B2 supergiants 

Prior to 2008, comparatively few spectroscopic data were avail- 
able for individual 09-B2 supergiants (Table 4), restricting anal- 
ysis primarily to the detection of line profile variability (LPV) 
in the Ha line. Of the stars with multiple spectra covering this 
transition, 12 have two or more observations for which R>2600, 
while a further 6 have two epochs of data of which one is of low 
resolution (R~500) and S/N. Of these 18, 9 stars demonstrated 
significant Ha LPV. This is illustated for W2a, 6a and 61a in 
Figs. 2-20, while it is apparent for W19, 23a, 28, 30a and 43a by 
comparison of the spectra presented in C05, Clark et al. (2008) 
and Negueruela et al. (submitted); the spectra of the final exam- 



5 Note that throughout the paper the rest wavelengths of transitions 
are indicated in Figures, while the stellar spectra demonstrate a systemic 
blueshift of ~40kms _1 . 
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Fig. 3. LPV in Ha for W61a (B0.5 la) and W70 (B3Ia). Note 
the variability in the strength of the temperature dependent C n 
photospheric lines in W70. 
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Fig. 5. Plot of I band magnitude versus spectral type for the late 
O - mid B stars within Wdl (stars individually dereddened ac- 
cording to the prescription given in Negueruela et al. subm itted] ! . 
The photometric variables identified by Bonanos (2007) are in- 
dicated by the filled symbols. 



pie - W28 - are omitted due to their comparatively low S/N. Poor 
temporal sampling places few constraints on the timescales of 
the variability, but in the case of W61a significant changes occur 
over a few tens of days (Fig. 3). Morel et al. (2004) monitored 12 
09-B2 Ia-Ib field supergiants, finding conspicuous variability on 
comparably rapid timescale in many of the stars with Ha profiles 
dominated by wind emission, including - as found for W19 - a 
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Fig. 6. LPV in Pal 1 indicative of photospheric pulsations for the 
BHGs W57a (B4 la) and W71 (B2.5 la). Similar variability is 
also present in the Ni lines. 



Fig. 8. LPV in the Ha line of W7 (B5 Ia + ) and W42a (B9 Ia + ). 
The insets show the Halpha profile with a log flux scale to high- 
light the broad emission wings. 
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Fig. 7. Plot of RV variations versus time for the Pal 1 line in the Fig. 9. LPV in the Ha for W13 (WNVL/B Ia + +OB) and W33 
BHGs W57a (B4 la) and W71 (B2.5 la). The RV is the error- (B5 Ia + ). 
weighted average of the Pa-11..14 lines, with an uncertainty of 
+3kms~ I . 



transition from absorption to full emission in HD 37128 (B0 la). 
Moreover all of the Ha line profile morphologies demonstrated 
by the OB supergiant population within Wdl were found to have 
identical counterparts in the samples presented in Kaufer et al. 
(fT996T > and Morel et al. (120041 . 

As such the variability of the Wdl members does not ap- 
pear atypical for late O/early B supergiants, while the poor sam- 
pling frequency of the current dataset suggests that many more 
stars are likely to be intrinsically variable. Such a supposition is 
supported by all 9 stars displaying Ha LPV being drawn from 
the subset of stars with 2 or more epochs of high resolution 
data. Although is not expected that such stars will have evolved 
far enough to encounter the LBV phase, the physical origin of 
such variability is currently uncertain, with pulsational, mag- 
netic and/or binary modulation of the wind all possible. In this 
respect it is of interest that W2a, 6a, 23a, 30a and 43a have all 
been identified as binaries (Clark et al. 120081 Negueruela et al. 
submitted and Ritchie et al. l2009al) . 




Fig. 10. Pulsational variability in the Paschen series and low ex- 
citation metallic lines of W7 over timescales as short as 35 days. 
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Of the 9 stars with both NTT (2002) and VLT (2004) I band 
spectra, only one star - W24 - showed evidence for variability; 
intriguingly this was most pronounced in the C m+Pal6 8500A 
blend rather than in other stronger, unblended Paschen lines. 
Motivated by this we examined the full VLT/FLAMES 2008-9 
spectroscopic dataset (with analysis of the comparable datasets 
for the full sample presented in a future dedicated paper). These 
data confirmed that the blend varied in both strength and radial 
velocity (RV) of the profile centroid over short timescales (~24d; 
Fig. 4). Given the weakness of the adjacent Pa 16 transition and 
Pal5+Can & Pal3+Can blends as well as the lack of compa- 
rable LPV in the Pal 1 line we associate this behaviour with the 
Cm 8500A line rather than Pal6 or Can,} 8498. The physical 
origin of this behaviour remains uncertain, with no periodicity or 
systematic secular evolution present in the current data. Ritchie 
et al. d2009al > report similar LPV for the Paschen series in the OB 
supergiants W8b, 21, 23a and 78, which they attribute to photo- 
spheric pulsations; it is tempting to include W24 in this category, 
although why pulsations should solely be visible in C m 8500A 
and not in the Paschen series is unclear. 

Another group of spectroscopically variable stars are the 
rare Of?p objects, in which the phenomenon may be linked to 
the presence of strong magnetic fields (e.g. Naze et al 2008). 
The defining feature of these stars is (variable) emission in the 
Cm/1/1 4647,50,51 multiple^ Unfortunately this feature is cur- 
rently inaccessible in W24, and while it is a hard X-ray source 
it is not over luminous for its bolometric luminosity as might be 
expected for an Of?p star (Clark et al. 2008). Nevertheless, given 
the presence of a magnetar within Wdl (Muno et al. 2006} and 
the suggestion that these are the descendents of highly magne- 
tised progenitors, the possibility that W24 could be a member 
of such a population is intriguing and warrants further investiga- 
tion. 

Regarding variability amongst the OB supergiant population 
within Wdl, we may also utilise the results of Negueruela et 
al. dsubmittedt to reinterpret the findings of Bonanos (20Q7J. In 
Fig. 5 we plot a semi-empirical HR diagram for all the OB super- 
giants within Wdl for which a spectral classification is available, 
highlighting the location of the photometric variables identified 
in the latter work. This emphasises that such stars are predom- 
inantly concentrated amongst the more luminous, evolved SGs; 
a similar result to that found by Fullerton et al. (1 19961 1. Given 
the apparent dependence on evolutionary state we assume that 
the photometric variability is most likely intrinsic to the star and 
hence not binary in origin. 

It is tempting to associate this behaviour with the a Cygni 
variables, of which the subset of late O/early B supergiants 
demonstrate quasi-periodic oscillations with <0.1mag amplitude 
(e.g. van Leeuwen et al. 1 1998b . Unfortunately, the observations 
of Bonanos (120071 1 - comprising two 10 day blocks separated 
by 20 days - are ill-suited to identifying the photometric modu- 
lation of a Cygni variables, which characteristically occur over 
2-6 week timescales for early OB supergiants, extending to hun- 
dreds of days for cooler super-/hypergiants (van Leeuwen et al. 
[19981 . 

Finally, four stars - W6a, 10, 1 1 and 28 - were also observed 
by West87; given their spectral classifications agree to <2 sub- 
types we find no evidence for their long term secular evolution. 
Nevertheless, of the 23 stars considered here, we find that all but 
four are either spectroscopic or photometric variables (Table 4). 



6 Note that the C ra 8500A line is the singlet counterpart (Ss'S-Sp'P ) 
to this multiplet. 



4. The mid-late B super-/hypergiants 

The 6 stars with spectral types B2.5-9 naturally divide into two 
subgroups. The first, containing W57a, 70 & 71, are earlier 
(B2.5-4 la) and have comparatively weak Ha emission, while 
the latter, comprising W7, 33 & 42a, are later (B5-9 Ia + ), and 
demonstrate stronger Ha emission with broad wings (likely due 
to electron scattering). Despite this, all were found to be spec- 
troscopically variable. 

Of the B2.5-4 la stars, the Ha profile of W70 was observed in 
absorption in 2002, full emission in 2003 and P Cygni emission 
in 2005 (Fig. 3). The C nvl/i6578, 6582 transitions also appeared 
to vary in both strength and the RV of the line centroid between 
observations, with changes in the profiles of both species ob- 
served over the shortest timescale probed (45 dfl In comparison, 
the early B supergiants HD 14134 (B3Ia) and HD 43384 (B3Iab) 
also showed a rapid (< 14 d) transition from Ha absorption to 
emission (Morel et al. 120041 ). although neither they nor the other 
B2-4 supergiants from that study demonstrated variability in the 
Cn lines. These transitions are temperature dependent and are 
absent in supergiants earlier than B0. 7 and later than B9 and 
peak at B3, although their apparent abundance dependance have 
precluded their use as an accurate diagnostic (Negueruela et al. 
submitted). Hence we may conclude that both wind and photo- 
sphere are variable, although we may not constrain the magni- 
tude of the temperature excursions exhibited by the latter. 

W71 also demonstrated significant LPV in Ha, which was 
present in emission in 2002 but absent in 2004, while it was 
absent in both spectra of W57a. For both objects the S/N and 
resolution of the spectra were too low to comment on the be- 
haviour of the nearby C n and He i features and for that reason 
are not shown here. Included in the binary survey (Ritchie et al. 
2009a), the multiple I band observations of both stars demon- 
strated rapid RV variability in the Paschen series and N i lines 
(Figs. 6-7). Building on these results we found ARV~12kms 
and 25kms~ 1 for W57a and 71 respectively with no period ap- 
parent in either dataset; consistent with random sampling of a 
short (pulsational) quasi-period, comparable in both respects to 
the photometric pulsations observed in other early and late B su- 
pergiants (Kaufer et al. 1997, 2006). The change in strength of 
the N i lines in both stars indicates that pulsations were present, 
being inconsistent with an origin for the RV changes in bina- 
rity alone. These imply moderate variations in the spectral type 
corresponding to AT~2000K, noting that modeling of the non- 
radial pulsations of the B supergiant star HD64760 revealed a 
comparable difference between minimum and maximum surface 
temperature (3500K; Kaufer et al. 120061) . 

The cool B hypergiants W7, 33 and 42a also demonstrated 
pronounced Ha LPV, although an emission component was al- 
ways observed (Figs. 8-9). The Cn doublet was also found to 
be variable in W7; this behaviour was mirrored in the Paschen 
series and low excitation metallic lines present in the I band 
spectra (Fig. 10). Variations in both the strength and RV of the 
C ii doublet have previously been observed for the B8 Iae star 
HD 199478 (Markovaet al. 120081) ; the similarity of the resultant 
spectroscopic period to the photometric period providing strong 
evidence for photospheric pulsations in this star. Identical LPV 
is also present in the I band spectra of W42a (not shown for 
brevity); in conjunction with comparable behaviour in both hot- 
ter and cooler supergiants within Wdl (W70 and 265 respec- 

7 While we regard the detection of RV changes as provisional, given 
the relatively low S/N and resolution of the observations, the changes 
in the line strength are robust, amounting to ~ 15% in Equivalent Width 
between the identical resolution 2005 FLAMES observations. 
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tively) we conclude that both W7 and W42a are also pulsating. 
As with W57a and 71, the range of wavelengths sampled by 
these spectra do not include particularly sensitive temperature 
diagnostics, although changes in spectral type by +1 subtype 
may be inferred for both W7 and W42 over timescales of ~35 
days (Fig. 10). 

However, as with the B2.5-4 la supergiants, the physical ori- 
gin of the Ha variability is unclear. Comparable P Cygni profiles 
and LPV in Ha w ere found for HD 199478, HD 92207 (AO la ; 
Kaufer et al. [19961 and the LBV HD 160529 (Stahl et al. [2003b . 
Both HD 199478 and HD 92207 undergo short lived (20-60 d) 
episodes of enhanced blueshifted absorption leading to the ap- 
pearance of a P Cygni profile (Kaufer et al. 119961 Markova et 
al. 120081 ): a timescale comparable to that inferred for the LPV 
in our stars (<45 days). Such 'high velocity absorption' (HVA) 
events are assumed to be due to a transient increase in the optical 
depth along the line of sight in an aspherical wind, although the 
origin of the asphericity is unspecified. However, in W7, 33 and 
42a the increase in strength of the P Cygni absorption trough was 
accompanied by an increase in the strength of the emission com- 
ponent of the line. Such correlated behaviour is not observed in 
HD 199478 and H D 92207, but is seen in the LBV HD 160529 
(Stahl et al. 120031) . where it is instead attributed to a long term 
secular change in the global properties of the stellar wind such 
as mass loss rate and terminal velocity. Further observations are 
clearly required to distinguish between both possibilities for the 
B hypergiants within Wdl 

In summary, evidence for both photospheric (pulsational) 
and wind variability are found for both sub-groups of B super- 
/hypergiants, although no causal link between both phenomena 
may be drawn from the current data nor physical cause for wind 
variability inferred. However, with regard to this we note that no 
binary markers are found for any of these stars; likewise there is 
no unambiguous evidence for either long term secular evolution 
or canonical LBV excursions over the 5-9yr observational base- 
line for individual stars (Table 5fl Such a result is not unprece- 
dented; P Cygni has demonstrated only weak secular evolution 
since its outburst in 17th Century - which would have been un- 
detectable in our data - with strong (>1 mag.) LBV excursions 
also absent (Lamers & de Groot 1992 de Groot et al. 12001b . 
Likewise the B hypergiants VI Cyg 12 and Sco - which show 
a remarkable spectral similarity to the BHGs within Wdl - have 
also been photometrically stable for a century (Clark et al. 2010, 
Sterken et al. [T99"7j l 

In relation to this, the spectral types determined by West87 
and C05 for W7, 42a, 70 & 7 1 are discrepant (Table 2). However 
given the difficulties of I band classification (Negueruela et al. 
submitted) we strongly suspect these arise from the low S/N and 
resolution of the spectra presented by West87 rather than from 
any intrinsic long term variability. In conjunction with photo- 
metric classifications (Table 2) we may therefore also exclude 
any long term transition to a significantly higher or lower tem- 
perature for these stars, such as the decades long excursions seen 
for R 127 and M33 Var A (Walborn et al. 120051 Humphreys et al. 
120061 Sect. 6). 

The final OB supergiant star we turn to is the sgB[e] star W9. 
Given their co-location in the HR diagram, there has been much 

8 I band spectra for LBVs in both high and low temperature phases 
are provided by Ritchie et al. (2009b ) and Munari et al. (f2p09 ), while 
observations of both early and late B/early A supergiants - the range 
over which a putative LBV within Wdl might be expected to evolve 
- are presented in C05 and Munari & Tomasella (1999 ). These data 
clearly demonstrate that changes in the spectral morphology resulting 
from a canonical LBV excursion would be identifiable in our data. 
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Fig. 11. LPV in Ha and Hei transitions for the candidate WN9 
binary W44. The inset in the left panel shows the same wave- 
length range but rescaled in flux to show the whole Ha profile. 



speculation on the evolutionary links between sgB[e] stars and 
LBVs. Both photometric (e.g. Bergner et al. 119951 Zickgraf et 
al. 11996) and spectroscopic (e.g. Zickgraf 2003 ) variability have 
been observed in these stars but to date these have not been un- 
ambiguously associated with characteristic LBV excursions or 
eruptions. In this respect W9 is of interest, given the evidence 
for changes in its unexpectedly high mass loss rate (Dougherty 
et al. 120101 ) and the detections of short term aperiodic photo- 
metric fluctuations (Bonanos 120071 ). A full description of the 
optical and IR spectroscopic dataset is presented in Ritchie et 
al. (submitted) but, perhaps unexpectedly, we found no evidence 
for pronounced spectroscopic variability over the course of the 
observations (2001-4; noting the comparatively low resolution 
data would not resolve LPV similar to that observed in other hot 
supergiants). 

5. The Wolf Rayets 

WRs have long been implicated in the LBV phenomenon; 
in their hot, quiescent state LBVs such as AG Carinae dis- 
play a WN9-11 spectral morphology, while the WN3 star 
HD 5980 underwent a major eruption resulting in a transition 
to significantly cooler spectral types in the 1980s (WN6 and 
B1.5Ia + ; Koenigsberger et al. 11998k Comparison of the optical 
and near-IR spectroscopic datasets (Clark & Negueruela 2002 
Negueruela & Clark 120051 Crowther et al. 120061 ) indicate that no 
such extreme behaviour is evident for any of the WRs within 
Wdl between 2001-5, with the discrepancies between the sub- 
types of individual objects being attributable to a combination of 
low S/N and resolution spectra and a poorly callibrated I band 
classification scheme. Unfortunately, the faintness of the WR 
population limits the utility of the long term photometry; of the 
8 stars observed on two or more epochal only W44 (=WR L) is 
potentially variable. 

As with OB supergiants, line profile variability has been ob- 
served in a large number of both WN and WC WRs (e.g. Robert 
119941 Lepine & Moffat [T999l ). typically taking the form of a su- 
perposition of variable subpeaks on an underlying broad emis- 
sion line and attributed to both small (stochastic clumping) and 

9 WRs A, B, E, F, G,L, M&S 
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large scale (co-rotating) wind structure. The origin of the latter 
is uncertain, with both photospheric and binary origins advanced 
for the quasi-periodic line profile variability observed in e.g. EZ 
CMa (=WR6; Morel et al. 120051 Georgiev et al. [19991) . although 
apparently unambiguous binary modulation has been reported in 
a number of WR colliding wind systems by Stevens & Howarth 
dT999T >. 

Including the extreme B hypergiants W5 (=WR S) and W13, 
which show a close continuity in physical properties to the WN9 
star W44, a total of five WRs have multiepoch optical spectro- 
scopic data, the others are the WC9 stars W239 (=WR F) and 
W241 (=WR E). Negueruela et al. (submitted) present the high- 
est resolution and S/N I band spectrum of W5 available to us; no 
evidence for long term line profile variability between 2001-4 
was visible in these data. W13 is discussed in depth in Ritchie 
et al. d2009aK who showed it to be a B Ia + /WNVL+OB binary. 
In addition to reflex RV motion in the Paschen series, significant 
variability was present in the emission components of profiles 
obtained at comparable orbital phases, behaviour replicated in 
the Ha profile (Fig. 9). Upon consideration of the full dataset no 
secular evolution was apparent in any spectral feature. 

Likewise, comparison of optical data from 2001-5 showed 
no long term changes for the WN9 star W44. However high re- 
sultion spectra revealed dramatic LPV in Ha and Hei AA 6678, 
7065 (Fig. 11) which, to the best of our knowledge, are more pro- 
nounced than the variations observed for any other WR star out- 
side of an LBV excursion. Given their absence in the spectra of 
our cool B hypergiants, the continued presence of strong He lA, A 
6678, 7065 emission throughout this period argues against at- 
tributing the changes to a transition to a cool LBV phase. On 
the basis of its X-ray properties W44 is a strong CWB candidate 
(Clark et al. 2008) and despite the above concerns - and lack of 
photometric confirmation - as with W13 we consider binary in- 
duced wind structure to be the most likely explanantion for the 
spectral behaviour, but further data are required to confirm this 
hypothesis. 

Finally we turn to the WC stars W239 and W241 (Fig. 12). 
Ritchie et al. (2009a) showed that W239 demonstrated reflex 
binary motion between 2008-9 but, while intrinsically variable, 
the expanded dataset revealed no evidence for secular evolution 
between 2001-9. Spectra of the dusty WC9 star W241 in 2005 
May and July demonstrated notable changes in the strengths of 
the Hei+n and Cii+iii lines and following the arguements for 
W44 we suspect these too are due to binarity. Of interest are a 
number of sharp emission features superimposed on the broad 
Cn emission line at ~6580A. We suspect these reflect emission 
in Ha and CnAA 6578, 6582, with a further, tentative identi- 
fication of emission in the [Nn]/1 6548 line also possible (the 
stronger [Nn]/16584 line is potentially blended with C nA 6582). 
The origin of such lines is currently unclear. We would not ex- 
pect narrow Ha or [N n] emission to be associated with a WC 
star, suggesting they must originate in the hitherto unseen com- 
panion. However, given the current data we are unable to deter- 
mine whether the emission lines are nebular in origin or instead 
whether they arise in a stellar wind; if the latter is correct they 
would imply a rather early spectral type for the secondary. 

6. The Yellow Hypergiants 

As with the hotter B hypergiants, field YHGs are also highly 
variable over multiple timescales. Most, if not all, demonstrate 
short period pulsations with periods of the order of 10 2 days (e.g. 
de Jager 1998 and refs. therein) visible in both lightcurve and RV 
variations in photospheric lines; similar rapid variability is also 
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Fig. 12. Variability in the spectra of the WC9d stars W239 and 
W241. Note the sharp emission features superimposed on the 
broad C n emission line of W241, with a possible identification 
of the [N u]A 6548 line in the blue wing (indicated in the inset). 



present in wind contaminated emission lines such as Ha (e.g. 
Lobel et al. 120031 ). Additionally, long term secular temperature 
variations have been inferred from both spectroscopic and pho- 
tometric observations for a number of the most luminous YHGs. 
Oudmaijer et al. (1996) found that IRC +10 420 increased in 
temperature by > 1000 K over the past twenty years, with similar 
behaviour inferred for HD 179821 (Patel et al. 120081) . In contrast 
HR 8752 has undergone two red loops in the HR diagram over 
the past thirty years, each of approximately 10 years duration 
(de Jager & Nieuwenhuijzen 1997). Similar excursions to cool 
temperatures have been observed for p Cas (Lobel et al. 2003 ) 
although the duration of these events is significantly shorter (<2 
yrs). 

Two or more observations of the Ha profile are available for 
5 of the 6 YHGs within Wdl. Of these we infer Ha to be vari- 
able in W265 (reported to be present by West87 but absent in our 
data), observe significant evolution in the profile of W16a (Fig. 
13; apparently absent in West87) and find no variability for W4, 
8a and 12a3 While only 3 epochs of spectroscopy are available 
for W16a, the profile appears to show a number of unusual fea- 
tures. The first, a narrow blueshifted emission component super- 
imposed on a presumably photospheric absorption profile is also 
present in t he YH Gs HD 179821 (Patel et al. 120081 and p Cas 
(Lobel et al. 120031 ) where, like here, it is variable. A broad plinth 
of emission extends to some ±1000 kms -1 and appears constant 
in all three spectra, while a third emission component extending 
from ~-700kms _1 - +300kms~ 1 seems to be present in the 2002 
spectrum but absent in 2005; to the best of our knowledge this 
combination of features is not present in any other YHG. 

Pronounced emission wings are present in the mid-late B hy- 
pergiants W7 and 42a, as well as lower luminosity late B/early 
A field supergiants (Kaufer et al. 1996). Typically, such features 
are thought to arise from electron scattering of line photons but 
no strong emission is present in W16a, nor in several of the stars 
studied by Kaufer et al ( 1996). Moreover the high velocity emis- 
sion plinth is asymmetric in W 16a and is more pronounced in the 
red, although this might be due to blended emission in the coin- 



10 Only a single epoch of high resolution and S/N data is available 
for W 12a - as such we would not expect to detect the subtle Ha LPV 
observed for e.g. p Cas (Lobel et al. 12003b 
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cident [Nn]/16583 and/or CnAA 6578, 6582 lines (noting that 
the later appear unlikely given they are only seen in the extreme 
early B supergiants). This feature is also reminiscent of the He i 
AA 5876, 6678 lines in W9 (Clark et al. l2008l . which Ritchie et 
al. dsubmittedl l tentatively attribute to the presence of a wind col- 
lision zone. However, there is no evidence that W16a is a CWB; 
therefore we are currently unable to offer a convincing physical 
explanation for either the line morphology or variability of Ha 
inW16a. 

Preliminary observations of the neutral metallic lines of 
W265 by Ritchie et al. (2009a) suggested that it was pulsat- 
ing; the inclusion of a second season of observations allows us 
to observe the turnaround in the pulsational cycle, implying a 
quasi-period of order ~10 2 days (Figs. 14-15). Unlike the late 
B/early A supergiants studied by Kaufer et al. (1997), we find 
a depth dependence for the velocity fields of the photospheric 
lines, with the higher excitation lines - which form deeper in the 
atmosphere - showing less variability than the lower excitation 
transitions, although all lines vary in phase. Comparison of the 
spectra taken at the extremes of this cycle to the classification 
standards of Munari & Tomasella ( 1999]Fig. 16) show it to vary 
between spectral types F0 (7200K) and F5 Ia + (6570K), spend- 
ing most time in the cool phase. 

The magnitude of this change is directly comparable to that 
observed for p Cas in observations between 1993-5 (AT~750K; 
Lobel et al. |1998l ), although the 2000-1 outburst of p Cas, which 
resulted in a greatly enhanced mass loss rate, was accompa- 
nied by a more extreme evolution of spectral type than observed 
here (early M supergiant; Lobel et al. 2003). Moreover, (i) the 
timescale of quasi-pulsations (mean value ~300 days; de Jager 
et al. U998l and refs. therein), (ii) the average ARV~18kms _1 for 
low excitation (l-5eV) photospheric lines (Lobel et al. 119941 ), 
(Hi) the development of excess absorption in the blueshifted 
wings of both the Ca n and Fe i lines during the hot phase (due 
to enhanced mass loss leading to an increased optical depth 
in an expanding atmosphere; Lobel et al. 1998 Ritchie et al. 
2009a) and (iv) the significant velocity stratification in photo- 
spheric lines originating at different atmospheric depths (Lobel 
et al. 120031) in p Cas are all directly comparable to the same phe- 
nomena in W265. Thus both stars appear to demonstrate similar 
pulsationally driven mass loss, presumably via the same physical 
mechanism. 

Unfortunately, the spectral datasets available for W4, 8a, 12a 
and 16a are of insufficient resolution and S/N to readily identify 
futher pulsators via RV shifts in the photospheric lines, although 
variations in the strength of neutral metallic lines are present in 
W4, and hence by analogy to W243 and 265, provide evidence 
for pulsations. Given the relatively poor temporal sampling of 
these observations and the finding that W265 appears to spend 
most of its time in the cool phase, we suspect that the remaining 
stars may also turn out to be variable. Nevertheless these data 
do permit accurate spectral classificatiorQ we found W4 (=F3 
Ia + ), W8a (=F8 Ia + ), W12a (=F1 Ia + ) and W16a (=A5 Ia + ) to be 
systematically later by 2-3 subtypes than found by C05 utilising 
lower resolution data. Given classification of A-G supergiants is 
primarily accomplished via the relative strengths of weak metal- 



11 Interestingly, some evidence exists for the presence of abundance 
anomalies in these data. Ni absorption features are significantly en- 
hanced in W12a with respect to W265; given little temperature depen- 
dence is expected for these lines between Fl-5 Ia + we suspect this is 
the result of chemical processing. Likewise photospheric C i absorption 
lines are present in all these stars but are absent in the cool-phase LBV 
W243 (Ritchie et al. 2009b I, suggesting C depletion in the latter object. 
A full analysis of these data will be presented in a future paper. 
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Fig. 13. LPV in Ha in W 16a (A5 Ia + ) between 2002-5. Note 
the match between the strengths of the DIB in the 2003 & 2005 
spectra implies that the change in strength of the narrow emis- 
sion and absorption components is real. The location of C nAA 
6578,82 and [Nn]A 6583 are also indicated. 



lie lines we conclude that this discrepancy results from the dif- 
ferent quality spectra available rather than being of astrophysical 
origin. A similar conclusion may also be drawn for the system- 
atically later classifications of W4, 8a, 32 and 265 by West87 
(Table 2). The classifications for W12a appear noticeably dis- 
crepant; we suspect this is due to an erroneously early spectral 
type in West87; the strength of the Pal8+Oi, Pal6+Can and 
Pal5+Can blends are all directly comparable to W265 in his 
data, as we find in ours. 

Consequently, we find no spectroscopic evidence for long 
term secular evolution of the temperature of any of the YHGs 
over the quarter of a century between 1981 -2005 . Such a conclu- 
sion is bolstered by the photometry available; the -magnitude 
visual brightening found for IRC +10 420 between 1930-70 
(Humphreys et al. 1992) would have been detectable if sim- 
ilar had occured for any YHG between 1966-2006 (Sect. 2). 
However, given the poor temporal sampling we may not exclude 
shorter events such as the red loops and outbursts which HR 
8752 and p Cas have been observed to undergo (de Jager 1998 
Lobel et al. 120031) . although long lived excursions to cool tem- 
peratures, such as seen for M33 Var A - for which a transition 
from an F to M supergiant, accompanied by a 3 magnitude de- 
crease in the V band, persisted for ~45years (Humphreys et al. 
2006) - would have been detected. 



7. The Red Supergiants 

Finally we turn to the RSGs W20, 26, 75 & 237, supplementing 
our sparse optical dataset (Table 5) with the near-IR classifica- 
tions of Mengel & Tacconi-Garman (2007). West87 report Ha 
emission in W20 & 26 and its absence in W237, with no data 
for W75. We fail to detect Ha in W20, but weak emission is 
present in W237, although at a level that could have been missed 
in West87 (Fig. 17); by analogy with W26, we consider it likely 
to be a blend of Ha+[N n]. Variable line emission in the Balmer 
series (and low excitation species) has previously been observed 
in the RSG+OB binary VV Cephei stars (Cowley [1969] ). where 
it is thought to arise in the region of the RSG wind ionised by 
the hot companion. 
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Fig. 14. Pulsations in the photospheric Mg i and Pal 1 lines in 
W265 
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Fig. 16. Comparison of spectra of W265 taken in the hot and 
cool phase to the classification spectra of Munari & Tomasella 
(TT999b . 
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Fig. 15. Plot of the evolution of the RV of selected photo- 
spheric line in W265; included are the upper photospheric 
lines of Fei (multiplet 60, Xlow ~2-2 eV), Mgi (multiplet 7, 
Xhw ~4.3 eV), Ni (multiplet 1, Xlow -10.2 eV) and the Pall 
line (xiow -12.04 eV). 



In contrast W26 continues to demonstrate strong Ha emis- 
sion in our data (Fig. 18), as well as forbidden line emission 
from [Nn] and [S m]; given the strength of the former with re- 
spect to Ha we suspect this emission arises in nitrogen enriched 
ejecta rather than a stellar wind. No changes are apparent in the 
strength of the Ha+[Nn] emission in W26 between 2001-6, al- 
though the missmatched resolutions mean that we are insensi- 
tive to subtle LPV. To the best of our knowledge only one other 
RSG - WOH G64 - is known to show such a pronounced nebular 
emission line spectrum (Levesque et al. 2009; see Appendix A 
for more details). 

No line diagnostics for spectral type are present within the 
wavelength regions sampled by our spectra and so we are forced 
to rely on the strength of the TiO bandheads at ~7050, ~8250, 
~8440 & ~8860A, which increase in strength with decreasing 



temperature (e.g. Levesque et al. 2005 Zhu et al. 1 19991 Ramsey 
et al. 1 198 11 >. Unfortunately, the low resolution and S/N compli- 
cate an absolute classification for the spectra from 2001, but we 



note that the bandhead appears significantly weaker in W26 than 
W20 and W237, suggesting an earlier spectral type for the for- 
mer, while the strengths of these features in W20 & 237 clearly 
differ (Figs. 17 & 18). 

Higher resolution spectra of the ~8860A bandhead for W20, 
W26 (both 2002 June) and W237 (2004 June; identical to those 
from 2005 March-July) are presented in Fig. 19. These indicate 
that W26 is of unambiguously earlier spectral type than W20 
and W237, which in these data - unlike in the 2001 spectra - 
are indistinguishable from one another. We therefore may con- 
clude that the temperature of one or both of these stars has var- 
ied between the two observations. Based on the callibration of 
the ~8860A bandhead by Ramsey et al. ([1981) we may assign 
spectral types of M2 la for W26 and M5 la for W20 and 237, 
with an uncertainty of +1 subtype in absolute rather than rel- 
ative classification. An additional spectrum of W26 from 2006 
February (Fig. 18) suggests a significantly later spectral type at 
this time in comparison to previous spectra (no earlier than M4 
la from the strength of the -8440A bandhead; Zhu et al. |1999) . 
while our sole spectrum of W75 from the same run appears to be 
much earlier than any other RSG (no later than M0 la; Fig. 17). 

Given the different near-IR diagnostics employed by Mengel 
& Tacconi-Garman (2007) a straightforward comparison of their 
classifications to those derived from these data is difficult. 
However, they find W20 and W26 to be of comparable spec- 
tral type (M5 la) and later than W237 (M3 la), with their clas- 
sification of W26 being compatible with ours from the previous 
month. 

We therefore conclude that W26 appears to unambiguously 
demonstrate changes in spectral type (M2-5 la), with pro- 
nounced variations in the TiO bandhead diagnostics between 
2001-6. W237 also appears to be spectroscopically variable 
based its classification relative to W20 and W26, with a compar- 
ison of the photometry presented by C05 and Bonanos (2007) 
supporting this assertion (Sect.2; we cannot exclude the possi- 
bility that W20 is also variable). Finally, we simply note the dis- 
crepancy between optical (<M0 la) and near-IR (M4 la) classi- 
fications of W75 between February-March 2007. 

Both photometric and spectroscopic variability may be 
present in high luminosity RSGs, although the latter is thought to 
be particularly rare (e.g. Levesque 2009). The extreme RSG VX 
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Fig. 17. Low resolution spectra of the RSGs W20 & 237 from 
2001 June and a medium resolution spectrum of W75 from 2006 
February, with an inset centred on the Ha+[Nn] blend in W237. 
Note that no emission is present in this feature in W20, while the 
S/N is insufficient to comment for W75. 




Fig. 18. Spectra of W26 from 2001 June and 2006 February re- 
spectively; note that the resolution of the 2006 spectrum has 
been artificially degraded to better match that of the 2001 data. 
The TiO bandheads are significantly stronger in the latter spec- 
trum, indicative of a later spectral type. The inset shows the re- 
gion of the spectrum from 2006 centred on Ha at the original 
resolution. 



Sgr demonstrates 'semi-regular variations' between M5.5-9.8 la, 
accompanied by AV>6mag (AI~1.2mag)over a ~732 day quasi- 
period (Lockwood & Wing 1982), while S Per transited from 
M4-7 la between 1971-4 (Fawley et al. 119771 ) and has also un- 
dergone pronounced photometric variability (~6mag) over the 
past century (Chipps et al. [2004). Significant evolution in this 
behaviour is also seen; between 1889-1919 VY CMa underwent 
irregular variations with AV< lmag. and a mean quasi-period of 
600 d, after which the excursions increased in both magnitude 
(AV~3mag) and duration (~1 100 d; Robinson [T970j ). 

Given the inevitable constrains imposed by our limited 
datasets, we conclude that the variability of W26 and W237 is 
entirely consistent with the behaviour of these extreme RSGs, 
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Fig. 19. High resolution I band spectra of the RSGs W20, 26 
(2002 June) and 237 (2004 June) covering the diagnostic TiO 
bandhead at 8860A. 



implying that they too are unstable. Given the quasi-periodic na- 
ture of the lightcurves of VX Sgr, SPer and VY CMa, Lockwood 
& Wing (1982) liken this behaviour to a high luminosity ana- 
logue of the pulsations that are present in Mira variables. The 
presence of substantial circumstellar ejecta around three of the 
four RSGs within Wdl (Dougherty et al. 120101 ) which, in the 
case of W26, appears to be nitrogen enriched, supports the as- 
sertion of Massey et al. (2007 ) that this instability may be asso- 
ciated with significant mass loss. Finally, as with both the BHGs 
and YHGs, the long term photometry (Table 2) precludes any 
long term excursions analagous to that undergone by M33 Var 
A. 



8. Discussion and concluding remarks 

Despite the disparate nature of the photospheric and spectro- 
scopic dataset accumulated for Wdl over the past 50 years, it 
presents a valuable resource to investigate the properties of the 
massive stellar population and in particular to characterise the 
prevalence and nature of stellar variability. This is of particu- 
lar relevance given the unique census of evolved stars within 
Wdl, comprising examples of all post-MS stellar types thought 
to originate from stars with M; m - ft - a j, ~30-40M Q . Given the limited 
sampling of the combined dataset it was somewhat surprising 
to find that stars in all evolutionary states were variable, with 
this behaviour being attributable to both evolving wind struc- 
ture^) and stellar pulsations. Indeed these phenomena appear 
to be commonplace, with 26 of the OB super-/hypergiants and 
all but one of both the cool hypergiants and the Wolf Rayets 
included in this study being either photometrically and/or spec- 
troscopically variable. 

Wind variability was most evident in the Ha line of the OB 
super-/hypergiants, with the line transiting from absorption to 
emission in W70, and pronounced P Cygni profiles developing 
in several other stars. Similar behaviour has been observed in 
field stars although, as with this population, the physical mecha- 
nism^) driving the wind inhomogineities are poorly understood. 
In this respect the detection of both photospheric pulsations and 
wind variability in the mid B supergiants W70 & 71 and hyper- 
giants W7 & 42a is of considerable interest and all four stars 
would benefit from further observations with a higher cadence 
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to investigate potential causal links between the two phenomena 
(c.f. Markova et al. 2008 ). Moreover, these results again empha- 
sise the uncertainty inherent in determining mass loss rates for 
OB supergiants from a single epoch of Ha spectroscopic data. 

Significant wind structure is also inferred for the WRs in- 
vestigated. Given the identification of binary markers such as 
the presence of hot dust, hard X-ray emission and, in the case 
of W13 and 239, reflex RV variations, we attribute this to the 
effect of the companion via wind/wind interaction, noting that 
no secular evolution is present in any of the stars. This is en- 
couraging since it suggests that detailed modeling of these vari- 
ation may yield the underlying binary properties of the systems 
(e.g. Stevens & Howarth 1999 ) enabling a comparison to the less 
evolved stellar populations within Wdl, hence illuminating the 
effects of binary interaction driven by stellar evolution. 

Wind diagnostics are more difficult to identify and study in 
the cool hypergiants. Nevertheless variable high velocity absorp- 
tion in the blue wing of the metallic photospheric lines is present 
in W265, indicative of changes in the optical depth and hence 
physical properties of the wind. By analogy to p Cas, we at- 
tribute this to pulsationally driven mass loss. Variability in the 
Ha emission line profile is inferred for W12a and directly ob- 
served in W16a, although the morphology of this line is complex 
and hence the physical origin of this behaviour is consequently 
uncertain. 

Variations in both profile strength and the RV of the line 
centroid were found for the photospheric lines of both the hot 
and cool evolved stellar populations, the combination indica- 
tive of pulsations rather than binary motion (e.g. Ritchie et al. 
2009a). Amongst the OB supergiants within Wdl, this behavour 
had previously been reported for BO. 5-2 la stars; our results ex- 
tend this to encompass the full spectral range present (09-B4 la). 
Utilising the new spectroscopic callibrations of Negueruela et al. 
(submitted) we find that the presence of photometric variability 
(Bonanos 2007 ) appears to be a strong function of evolutionary 
phase, being more prevalent amongst the early-mid B la stars in 
comparison to the late O Iab-II stars. This is qualitatively similar 
to the conclusions drawn by Fullerton et al. ( 1996) regarding the 
occurence and magnitude of photospheric pulsations in OB stars 
and hence it is tempting to attribute the photometric variabil- 
ity to photospheric pulsations, such as those inferred for the a 
Cygni variables. Unfortunately, we currently lack the high S/N 
and cadence observations that would enable pulsational mode 
identification. 

Amongst the cooler super-/hypergiants, pulsations were in- 
ferred for stars with spectral types ranging from mid B through 
to mid F. The best studied example - W265 - was observed 
to transit between F0-5 Ia + over the course of a single cycle 
(~100 days); the combination of both timescale, magnitude and 
atmospheric depth velocity stratification of the pulsations being 
strongly reminiscent of those observed in p Cas, which have 
been attributed to non-radial modes (Lobel et al. 119941 12003). 
Although the data are sparse, both W26 and 237 display spectro- 
scopic variability indicative of changes in spectral type, with the 
former being observed to evolve from M2-5 la. The finding that 
at least two of the four RSGs within Wdl are unstable and by 
analogy with other extreme RSGs, possibly pulsating is of con- 
siderable interest given (i) their extreme luminosity and hence 
high progenitor mas£3> (H) their rich circumstellar envelopes, 
which provide evidence for significant recent mass loss, possi- 



12 A lower limit to their progenitor masses may be inferred from the 
cluster Main Sequence turnoff. This is an important consideration given 
that the mass/luminosity degeneracy for RSGs as they execute a red 



bly driven by these instabilities and ( Hi) the rarity of this phe- 
nomenon (with only six other examples known in the Galaxy 
and Magellanic Clouds; Levesque 2009]). 

Indeed, these preliminary findings emphasise the important 
role studies of young massive clusters with multiplexing spectro- 
graphs will play in understanding the physics of high mass stars. 
In particular the detection of photospheric pulsations extending 
from spectral types as early as 09 I, through F5 Ia + hypergiants 
and potentially as late as the M5 la supergiants confirms the 
finding of e.g. Burki (1 19781 1 that stars in this region of the HR 
diagram are unstable, but critically for a rich co-eval popula- 
tion of stars at a single metallicity rather than a heterogeneous 
field population. Given the sample size - ~80 ~08III-B2 la stars 
- future analysis of the full 2008-9 VLT/FLAMES dataset will 
provide valuable observational constraints on the location of the 
high temperature boundary of this region of instability in a sta- 
tistically robust manner. 

Furthermore, observational determination of the amplitude 
and period of the pulsations as a function of temperature, lumi- 
nosity and evolutionary phase will provide direct contraints on 
the underlying physical mechanism(s) driving this instability - 
for example strange- and/or g-modes (Kiriakidis et al. 1993 Saio 
et al. 120061 1 - and, in conjunction with chemical abundance anal- 
ysis, a detailed picture of the evolution of the internal structure 
of massive stars as they execute a red loop on their way to be- 
coming Wolf Rayets. Similarly, high candence observations will 
enable an investigation into the role pulsations play in both the 
origin of wind structure in hot supergiants and in driving mass 
loss in the cool hypergiants. Indeed, the uniquely rich population 
of both hot and cool super-/hypergiants within Wdl holds open 
the possibility of investigating the varying contribution of both 
line and pulsational driven mass loss as a function of stellar tem- 
perature and luminosity as stars evolve across the HR diagram. 

Relating to this, with the exception of the known LBV W243 
(Ritchie et al. 2009b), we find no compelling evidence for other 
LBV like excursions or eruptions within Wdl, despite both the 
large number of transitional stars present and the evidence for 
past episodes of enhanced mass loss amongst a subset of both 
hot and cool stars (Dougherty et al. 1201 Oi l. Ha LPV similar to 
that observed in the LBV HD160529 was found for several of 
the B super-/hypergiants, but the poor temporal coverage meant 
that we could not determine whether these corresponded to the 
long term secular changes that characterise LBV excursions, or 
the short lived variations indicative of an asymmetric wind. 

However, while similar limitations also apply to the (non- 
) detection of events comparable to the outbursts of p Cas and 
the red loops of HR8752 in the YHG population, we would 
have expected to detect events comparable to the multi-decade 
excursions of e.g. M33 Var A, or the long term secular evolu- 
tion of IRC +10 420. Given these results, it is of interest that 
a comparable analysis for the WNLh stars - which have been 
proposed to be quiescent LBVs - within the Galactic Centre and 
Quintuplet cluster similarly reveals little evidence for significant 
long term variability (Appendix B). Apparently transitional stars 
co-located with both 'bona fide' LBVs and rapidly evolving cool 
hypergiants in the HR diagram (e.g. Clark et al. 2005b) can en- 
joy significant periods of quiescence. 

These results are surprising given that the majority of the 
stars observed within Wdl, from the extreme-B hypergiants 
through to the RSGs, appear to be low-level variables and only 
borderline-stable. As such one might naively have expected them 



loop prevents such a determination for field stars (Meynet & Maeder 

HESJ. 
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to be prey to such large scale instabilities, which presumably 
contribute to the lack of stars above the HD limit. In this respect 
the YHGs and RSGs are of particular interest since they are lo- 
cated at the low temperature boundary of the HD limit and hence 
appear to represent the most massive stars that may exist in this 
region of the HR diagram at solar metallicities. 

Drout et al. (120091) present detailed evolutionary predictions 
for the physical properties of such stars, which we may test 
against the population within Wdl. Our new spectral classifica- 
tions confirm the expectation that YHGs evolving from ~4OM 
stars should have relatively high temperatures (> 6000 K), while 
the lack of secular variability for any of these stars over a period 
of ~50yrs is consistent with the expectation of a relatively long 
lifetime for this phase (~10 4 -10 5 yr). However at these luminosi- 
ties YHGs are not expected to evolve to cooler temperatures, but 
exactly this evolution is indicated by the presence of 4 RSGs 
within Wdl. An extention of the redwards loop to encompass a 
RSG phase would result in a corresponding reduction in the time 
spent as a YHG While rapid evolution through the YHG phase 
might be consistent with the real time increase in stellar temper- 
ature observed for IRC +10 420 and HD 179821, the continuing 
absence of this behaviour from the cool hypergiant population 
of Wdl will place increasingly stringent constraints on such a 
scenario. 

Therefore a synthesis of ( i) the current physical properties of 
the evolved stars within Wdl - such as chemical composition, 
luminosity, temperature and, where observable, surface gravity, 
( ii) the relative populations of the different evolutionary classes 
and (Hi) the characterisation of their (pulsational) instabilities 
and secular variability has great potential for constraining the 
physics of massive post-MS evolution. 

Finally, regarding the ensemble properties of Wdl, con- 
straining the occurence and amplitude of the pulsationally (and 
binary) driven RV changes will be crucial in allowing the ac- 
curate determination of dynamical masses for unresolved extra- 
galactic clusters (e.g. Gieles et al. 2009), where the velocity dis- 
persion attributable to such processes may match or even domi- 
nate virial motions. Likewise determining the properties of pho- 
tometric outbursts driven by pulsational instabilities in both hot 
and cool hypergiants will be important in determining the va- 
lidity of photometric cluster mass estimates, particularly for the 
lower range of masses where the presence of such an event may 
dominate the integrated light. 
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Appendix A: The emission line spectrum of the 
RSG W26 

Of the four RSGs within Wdl, only W26 demonstrates a pro- 
nounced nebular emission line spectrum, with emission in Ha, 
[Nn] ,l/i6548,6582 and [Sm] ^9069,9532. The [Sra] lines are 
common in high temperature Hn regions as well as tracing 
shocked material. No oxygen emission - such as [Oi] A 6300 
- is seen, while [S in] A 6312 is particularly temperature depen- 
dent (e.g. Osterbrock & Ferland 2006 ) and so its absence while 
[S in] /Li9069,9532 emission is present is explicable if the re- 
gion is not highly excited. The only other RSG to demonstrate a 
nebular emission line spectrum is the extreme system WOH G64 
which, like W26, is also associated with an extensive dusty cir- 
cumstellar envelope (Levesque et al. [2009). Unlike W26, WOH 
G64 lacks [S m] emission although [Oi] A 6300 and [S n] A 6713 
are present, as is the high excitation line [O m] A 5007, which 
our spectra do not reach (Levesque et al. 20091. 

Given the lack of many (blue) nebular diagnostic lines we 
are unable to present quantitative conclusions as to the chemi- 
cal abundances or excitation mechanism for the emission spec- 
trum. However, by analogy to both WOH G64 and the nebula 
around the BSG Sher25, we suspect the strength of the [Nn] 
emission implies N enrichment in the circumstellar material, 
presumably as the result of CNO processing (Levesque et al. 
2009). Likewise, while we suspect the [S m] emission is shock 
excited, we cannot exclude the possibility of excitation by either 
a hot companion or the diffuse UV radiation field of the cluster; 
indeed the location of W26 within the 'core' of Wdl suggests 
that the latter may well play a significant role (c.f. Dougherty et 
al. 120101 ). Nevertheless, the lack of comparable nebular emission 
in W20 and 237 is notable, particularly given that these 2 RSGs 
also possess extended circumstellar envelopes comprising both 
dust and ionised gas; albeit significant less massive and dense 
than found for W26 (Dougherty et al. l20T0l ). 



Appendix B: Long term stellar variability in the 
Galactic Centre and Quintuplet clusters 

Comparable variability analyses using published data may also 
be undertaken for both the Galactic Centre (GC - 6 Myr; 
Paumard et al. 120061 ). and Quintuplet clusters (~4Myr; Figer 
et al. I1999I ). Multiple spectroscopic observations spanning a 
decade have been made of the emission line star population 
within the GC cluster (e.g. Najarro et al. 119971 Paumard et al. 
120011 Martins et al. 2007 ) which are complemented by the long 
term (~decade) photometric datasets of Rafelski et al. (2007) 
and Trippe et al. ((2006 ) . In the absence of cool hypergiants and 
comparable multiple spectroscopic datasets for the OB super- 
giant population we restrict the discussion to the narrow line 
He i stars, which recent analysis by Martins et al. (2007!) con- 
firm to be Ofpe/WNL stars. Despite the implication of such 
stars in the LBV phenomena only one star - IRS34W - was ob- 
served to undergo photometic variations of a magnitude com- 
parable to canonical LBV excursions (Trippe et al. 2006 ), with 
a second, IRS16NW, being a low level variable (AK~0.2mag). 
However, with the exception of reflex radial velocity motion in 
the eclipsing binary IRS16SW, no spectroscopic variability was 
found for any of these stars, including IRS34W (Trippe et al. 
2006 1 Martins, priv. comm. 2009), noting that the S/N and reso- 
lution of the spectra would clearly have identified any significant 
changes in morphology characteristic of the LBV phenomenon 
(e.g. Clark et al. 120091 ), as well as the LPV exhibited by e.g. W7 
and W33. 

Regarding the Quintuplet, both the cool B hypergiants 
FMM362 and the Pistol star are spectroscopically (Figer et al. 
1999; Figer, priv. comm. 2009) and photometrically variable 
(Glass et al. 1999), suggesting that they are undergoing LBV 
excursions. Unfortunately, only two epochs of spectroscopy are 
available for the remaining cluster population, separated by of 
order a decade (Figer et al. 119991 Liermann et al. 2009). Of these, 
the former are of significantly lower resolution and S/N, al- 
though are still sufficient to distinguish between e.g. cool hyper- 
giants such as the Pistol star, hotter Ofpe/WNL and ~featureless 
OB supergiant spectra. However, with the exception of FMM362 
and the Pistol star, none of the other stars in common between 
the studies appear to demonstrate unambiguous spectral vari- 
ability, despite the presence of a number of luminous WN9 and 
late O/early B supergiants within the Quintuplet. These results 
are bolstered by Glass et al. (1999) who only found two other 
confirmed cluster members to be photometrically variable be- 
tween 1994-7; the dusty WCL stars Q2 and Q3, for which this 
behaviour may be attributed to time variable dust production. 

Thus, somewhat suprisingly, of the 10 WN9 and WN9/Ofpe 
stars for which multiepoch data exist, none have demonstrated 
characteristic LBV spectral variability, while only a single object 
would have been photometrically classified as such. Mindful of 
the sparse nature of the Quintuplet datasets, these findings are 
suggestive of a rather long duty cycle for such excursions if such 
stars are quiescent hot phase LBVs. 



